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Abstract: We present a comprehensive study of the cosmological solutions of 6D 
braneworld models with azimuthal symmetry in the extra dimensions, moduli stabiliza- 
tion by flux or a bulk scalar field, and which contain at least one 3-brane that could be 
identified with our world. We emphasize an unusual property of these models: their ex- 
pansion rate depends on the 3-brane tension either not at all, or in a nonstandard way, 
at odds with the naive expected dimensional reduction of these systems to 4D general 
relativity at low energies. Unlike other braneworld attempts to find a self-tuning solution 
to the cosmological constant problem, the apparent failure of decoupling in these models 
is not associated with the presence of unstabilized moduli; rather it is due to automatic 
cancellation of the brane tension by the curvature induced by the brane. This provides 
some corroboration for the hope that these models provide a distinctive step toward un- 
derstanding the smallness of the observed cosmological constant. However, we point out 
some challenges for obtaining realistic cosmology within this framework. 
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1. Introduction 

The cosmology of codimension one branes, notably 3-branes in a 5D universe, has been 
intensively studied and is now rather well understood. One of the striking predictions of 
this picture is that the Friedmann equation for the Hubble expansion rate should have the 
form = ^Gp{l + ^) in the simplest situation with a single brane of tension T and 
excess energy density p, in a warped space where the bulk is 5D anti-deSitter space Q-lll. 
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One naturally wonders whether this kind of behavior is particular to branes of codimension 
one {i.e., having only one transverse dimension), or if higher codimension branes can act 
similarly. There have been numerous attempts over the last twenty years to harness some 
of the unusual features of six dimensional models in order to get some insight into the 
cosmological constant problem |Q]-[14|, the weak scale hierarchy problem |5|, |6[,|jl5|-[^ or 



other effects of warped compactification [23|-|28], and moduli stabilization p^,|29|-p2|. 
However we are not aware of any attempts to systematically explore modifications to the 
Friedmann equation in codimension two braneworld models. In fact, most of the work 
referred to above has been confined to the study of static solutions. 

In this paper we consider the most symmetric case, where the extra dimensions have 
azimuthal symmetry and the metric functions depend only upon the radial coordinate of 
the bulk, with periodic angular coordinate 9 G [0, 27r]: 

ds"^ = M^{p)g^^dx^dx'' + dp^ + L'^{p)de^- (1.1) 

here g^u is a maximally symmetric 4D metric which can be deSitter, anti-deSitter, or 
Minkowski space, with cosmological constant A4, normalized such that the scale factor is 
a{t) = e^* with H = KgyAl/S in terms of the 6D gravitational constant kq. General FRW 
solutions in the 4D spacetime would require having distinct warp factors for dt'^ and (dx*)^, 
a complication which we choose to avoid in the present study. We will assume that the 
only bulk contributions to the stress energy are a 6D cosmological constant Ag, a magnetic 
flux Fpg, or a bulk scalar field (j). 

In section 2 we will classify all the possible bulk solutions which are consistent with 
( |1.1| ) and point out the insensitivity of A4 to 3-brane tensions in some of these solutions, 
including a particularly novel one in section 2.1. In section 3 we discuss the boundary 
conditions associated with the tensions of branes which bound these solutions, which are 
the source of any dependence of the Hubble rate on the brane tensions which does exist. 
Section 4 sets the stage for investigating the Friedmann equation in a class of warped 
solutions by reviewing how these solutions can be stabilized by a bulk scalar field. This 
solution is known as the AdS soliton, since it looks like AdSg at large radii. In section 5 we 
derive a surprising result for the AdS soliton, namely that the Friedmann equation relating 
the expansion rate to the brane tensions has an unusual form which does not agree with 
the expectation that the system should reduce to 4D general relativity at distances much 
larger than the compactification scale. In section 6 we turn our attention to an unwarped 
solution, studied recently in ||l^, U], where the bulk is a two-sphere stabilized by magnetic 
flux and deformed by deficit angles due to 3-branes at antipodal points. We look for 
unstable modes in the spectrum of fluctuations around this solution which might explain 
why the expansion rate of these models is insensitive to the 3-brane tensions, and show that 
no such modes exist. We conclude with a discussion of how 2D compactification manifolds 
manage to tune away the effects of 3-brane tensions in certain cases, and the difficulties 
that will have to be overcome if one wants to study realistic cosmology by putting matter 
with a general equation of state on the branes. 



- 2 - 



2. Bulk Solutions 



We start by writing the Einstein and scalar field equations for the metric ( |1.1D . Defining 
fi = M' /M, i = L' /L, and working in units where the 6D gravitational constant Kq = 1, 
these equations are 

M^: ^' + 3^' + ^2 + 6M^ + 3£/. = -A6 + ^-^(<A'^ + mV + ^)-^ (2.1) 

9 . 2A4 1 / . , n^(p'^\ , , 

■ V + 10/^^ = + - 2 + ^ - + 21^8 (2-2) 

2/2 

0: + (4^ + = m2</> + ^ (2.4) 

The 4D cosmological constant is defined by = A4/3 in the de Sitter brane case, where 
the 4D line element has the FRW form ds^ = —dt^ + e^^^dx'^ We have allowed for the 
scalar field to be either real or complex with a winding number n. The constant /3 is related 
to the magnetic field strength by = PL/M'^, which can be seen to satisfy the Maxwell 
equation dA{\/\G\F^^) = 0. For simplicity we have assumed that the bulk scalar has no 
self interactions, so V{(l)) = ^rn^cf)^. Only three of these equations are independent; for 
example the (/i^) equation can be derived from differentiating {pp) and combining with 
the other equations. 

Rubakov and Shaposhnikov found an elegant way to solve the 6D Einstein equations 
resulting from the metric ( |1.1D , which was generalized to include bulk magnetic fiux in Q 
(see also [||). We extend the method to take into account the possible presence of bulk 
scalar fields. The Einstein equation (^]^) for this system is equivalent to the equation 
governing the classical motion of a particle of unit mass, whose action is 

S = j dp - U{z) - U^{z)^ . (2.5) 

Here z' = ^ and the potential depends on the 4D and 6D cosmological constants and the 
magnetic field via^ 

C/(z) = az2-?,z6/5^^^-6/5. a = ^k,- 6=^A4; c = ^/J^ (2.6) 

lb 24 9b 

[/,(.) = |z2(<A'^ + mV-^) (2.7) 

The radial coordinate p plays the role of time in this analogy. If we ignore the bulk scalar 
(j), the resulting equation of motion z" + U'{z) = allows one to solve for the warp factor 
M{p) independently of any other fields in the problem, using the relation 

M{p) = z^l\p). (2.8) 



^we correct a factor of 2 error in 6 relative to 
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If (j) is nonzero, this is no longer possible; however one can find an approximate solution 
iteratively by first finding M{z) when i;^ = 0, next solving for (j) in the background geometry, 
and then finding the approximate form of the back reaction of (f) on the geometry by treating 
as a perturbation. 

The above procedure has so far yielded only the metric component M. To obtain L, 
we can in general solve the first order {pp) equation; but in most cases there is a much 
simpler relation which gives either an exact expression or a good approximation to L in 
terms of M: 

L(p) = r'-^. (2.9) 

dp 

Here i? is a constant of integration, which as we will show later is determined by the tension 
of a 3-brane which may be located at p = 0. This result follows from the difference of the 
{09) and {pp) equations, which can be written as 

L' M" , , 

— = + — (2.10) 

L M' Ap ^ ' 

Thus the relation ( |2.9D is exact whenever cf)' = with the exception of unwarped solutions 
where M' = 0. For these the {pp) Einstein equation provides no information about L, and 



we must turn to the {pp) equation, as will be discussed in section 2.3 . 

Let us now focus on situations in which (j)' is negligible and the scalar field has no 
winding; in this case can be absorbed into the a term of U . The potential U{z) has 
eight possible distinctive shapes depending on the signs of p and A, and whether or not 
there is a magnetic flux. These are summarized in figure 1. The solutions can be visualized 
by starting the particle at rest {z' = 0) at position p = with some initial value zq, and 
letting it roll in the potential. The position p should be thought of as time in the mechanics 
analogy. The initial condition that z' = implies that M' = at /> = 0, which must be the 
case; there is no singular source term for M" in the Einstein equations, and therefore M' 
must be continuous at the origin. Considering M' evaluated along a straight line passing 
through the origin, one sees that continuity of the slope requires that M' in fact vanishes. 



integrating the equation of energy conservation, E = iz'^ + U{z): 



Using the mechanics analogy, a closed form solution for M{p) is implicitly provided by 

1 . 

dz 



P 



^2{E - U{z)) 



(2.11) 



Once the particle starts rolling, there are four possible outcomes, which depend on the 
shape of the potential and the initial value zq. 

1. In one kind of solution, the particle can reach z = at some finite value of p. This 
happens when (3 = 0, Aq > and E > (fig. la), or when /3 = 0, Ag < 0, A4 < 
and z < Zm (fig- lb). Some curvature invariants diverge, so the extra dimensions end 



at a singularity with the topology of Si in this case [23]. 



2. Another kind is where z ^ 00 as p ^ 00, as occurs if Ag < and A4 > 0, or if 
Ag < 0, A4 < and z > z^. Then gravity will not be localized, and the world does 
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Figure 1: Shapes of the potential U{z) depending on the sign of A4 for (a) Ag > and (b) Ag < 
in the case of vanishing magnetic field (3, and similarly in (c) and (d) for nonzero f3. Numbered 
dots show initial conditions leading to the four kinds of solutions discussed below. The first critical 
value is given by A^^ = {2Ae/5)^P'^. 



not look 4-dimensional to an observer on the brane. A special case of this kind is the 
AdS soliton solution pJ, where /3 = A4 = 0. 



A third possibility is that the particle starts to the left of a stable equilibrium position, 
in which case it comes to rest again at the other side of the minimum at some 
maximum value of p = pm- The vanishing of z' means that the space closes off at pm, 
and we can consistently insert a second 3-brane at this position. The extra dimensions 
look like S2 deformed by deficit angles due to the two 3-branes at antipodal points. A 
special case 14] is when the particle stays at rest at the equilibrium position. This 



gives an unwarped solution since M(p) remains constant. Nonvanishing magnetic flux 
/3 is required to obtain static solutions of this kind. 

4. One can also find special unwarped solutions when the particle sits at an unstable 
equilibrium position, denoted in fig. 1 by Zm, in the case of an AdS 3-brane. Since 
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we are interested in cosmological solutions, we do not pursue these here. 

To cure the problems of cases (1) and (2) above, one should cut out the part of the 
space with z = or z = oo at some value z = Z4, and insert a 4-brane at the corresponding 



radius p = Pm [22|. In fact it is not generally possible to satisfy the jump conditions (to be 
described in section 3) at pm using an ordinary 4-brane, since the bulk solution requires 
that its stress-energy components Tgg differ from the other spatial components Ta. One 
way of accomplishing this is to imagine that a pure tension 4-brane is accompanied by a 



3-brane which is smeared around the compact extra dimension |34], or by Casimir energy 



of a massless field confined to the 4-brane [ 17 1 . Alternatively, it is possible to find a special 



value of pm where the jump conditions on M and L can both be satisfied for a pure tension 



4-brane, as long as A4 is nonzero. By integrating the exact result ( 4.10 ) which we will 
derive in section 4.3, this value is implicitly given by 



p™ . L'{0) 
A 



M^Ldp=^^ (2.12) 



showing that pm ^ 00 as A4 — > for the warped solutions where M and L grow expo- 
nentially. If one wants a solution which remains compact in the static limit, this is not 
acceptable. 

For completeness we mention one other kind of solution: we can start with z' ^ 
at p = 0, but in this case the position p = should not be regarded as a single point, 
but rather the location of another 4-brane. A special case is when Ag < 0, A = = 0; 
then the warp factor takes the simple exponential form of the RS model. In fact this is 
just the 5D RS solution augmented by one extra compact dimension, whose warp factor 
is identical to that of the large dimensions. This model does not involve 3-branes, and is 
mathematically quite similar to the 5D case which has been so well studied already. A 
more thorough examination of the cosmology of these solutions in presented in appendix 
B, and we comment upon its qualitative differences relative to the codimension two case 
in the conclusions. 

In the following subsections we present some specific analytic solutions illustrating the 
cases mentioned above. 

2.1 Type 1 solutions 

A solution which is borderline between type 1 and type 3 can be found by choosing A4, Ag > 
0, E = 0. Starting from z = Zc = (b/a)^^^, the critical value shown in fig. 1(a) defined by 
U{zc) = 0, the particle subsequently rolls toward z = 0. The solution is 



M{p) = cos (kp); L{p) = Rk sin {kp); k = \ = H (2.13) 




where we have used the freedom to rescale x^^ to make M(0) = 1. Recall that H = y^W^ 
is the Hubble constant, so that A4 is determined by 

A4 = ^Ae (2.14) 
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Interestingly, no curvature invariants diverge at the point pm = 7r/2fe where M = and 
the space should be terminated. In fact the Ricci scalar is constant, R = SAg, and so are 
the Ricci tensor squared RapR"^^ = |Ag, and Riemann tensor squared Rap^sR^^'^^ = f-'^e- 
The Weyl tensor vanishes. This is analogous to the behavior as y ^ oo in the 5D Randall- 
Sundrum solution [^] with the line element ds^ = e~'^^^y\—dt^ + dx"^) + dy"^; even though 
the warp factor vanishes as y ^ oo in this model, there is no singularity. Although 
M'{pm) / 0, yet z' = b/2M^/'^M' vanishes at Pm, so in the mechanics analogy, the ball 
rolls to the top of the hill and comes to rest again. Apparently the 4D part of the spacetime 
disappears at pm, leaving one with just 2D Euclidean space in the p-6 plane. 

In the case of vanishing A4 and positive Ag, the following solution can be found P] 



M{p) = cos{kpf/^ ■ L{p) = Rksin{kpf'^; k = yj^KQ. (2.15) 

In this case however, in the mechanics analogy the particle does not come to rest at z = 0, 
and there will be a curvature singularity at that point. It is therefore necessary to insert a 
4-brane &t p < ps to cut off the space before reaching the point where it becomes singular. 

2.2 Type 2 solutions 

It is interesting to consider solutions with Ag < since this gives a warped geometry, 
where the hierarchy problem can be solved on the 3-brane d la Randall-Sundrum. We find 
a solution of type 2 when P = E = and A4,A6 < by analytically continuing ( 2.13| ) in 
A4 and Ag: 



M(p) = cosh (kp) ; L{p) = Rk sinh (kp) ; A; = ^ ^ = iH (2.16) 

This follows the upper curve of fig. 1(b), starting at Zc- The Hubble constant is imaginary 
because the 4D metric is AdS space. In addition the E = type 2 solution with A4 > 
and Ag < is 

M{p) = sinh [kp) ; L{p) = Rk cosh [kp) 

which follows the lower curve of fig. 1(b), starting at z = 0. However, this solution does 
not admit a 3-brane at p = since L is nonvanishing at this point. 

The only way to admit a 3-brane at p = if A4 > is to take ii^ < 0, so that the 
initial condition z' = (equivalent to the boundary condition M' = 0) can be satisfied. 



The static solution with A4 = (and /? = 0) has been investigated in [p4| , 22|; this is the 
AdS soliton truncated at large /o by a 4-brane: 



M(p)=coshf(M; L{p) = -Rk^^^^ k = J-lK, (2.17) 

cosh5(A;p) V » 

Analytic expressions for M{p) in the AdS soliton type of solution no longer exist when 
A4 > 0. However, we can obtain approximate analytic solutions by considering A4 as a 
perturbation. The details are given in appendix A. The result is that the unperturbed 
solution gets shifted according to 

M{p) = [zo cosh{k{p + 5p))f^ ; L{p) = RM'{p) (2.18) 
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where 6p{z) has the shape shown in fig. 2, and z = zq cosh{kp). At sufficiently large values 
of p, the shift 6p approaches a constant, so at leading order the only difference between 
the perturbed and unperturbed solutions at large p is a shift in the radial size of the extra 
dimension. 




1 2 3 4 5 6 7 8 9 



z 



Figure 2: Sp versus z, due to nonzero A4 in the AdS soliton solution. 



2.3 Type 3 solutions 

The simplest example of a solution of type 3 is that in which the particle sits at a minimum 
of U. In this case M is constant, so the solution is unwarped. If there is no magnetic 
field, this occurs only in fig. 1(a), when A4 > 0. To specify the physical value of A4, we 
should rescale M ^ 1 so that time and energy are normalized in the usual way (ds^ = 
Since A4 = which has dimensions of (mass)^, the relevant quantity 
is K^/M'^ = A6/2; in other words the expansion rate is governed completely by Ag, which 
is to be expected since it is the only source of stress energy. 

Another example of unwarped solutions is the static case A4 = with nonvanishing 
magnetic field. Just like above we noted that the physically relevant combination for the 
expansion rate is Ai^jM'^ ^ in this case fp' jM"^ is the meaningful combination for the square 
of the magnetic field strength. Minimizing U{z) we find that 0^ /M^ = 2Ae, in agreement 
with references [13, 14 1. Notice that Ag must be positive to achieve this fine tuning with 
the magnetic flux. 

In fact, the general relation between A4, Ag and (3 can be easily found in the same 
way, demanding U\zq) = 0, when we realize that there is the freedom to set zq ^ 1 by 
rescaling x^/zq. Thus 

1 1 7 

(2.19) 



A4 = ^Ag - 



characterizes the physical expansion rate for an observer on the 3-brane when the bulk 
parameters are not tuned to give a static solution. This is a significant result because it is 
completely independent of the tension of such a 3-brane |l^, |l^ . 

Although M is trivial for the unwarped solutions, L is nontrivial: 

L{p) = Rsinikp); fc^ = ^Ag + ^/^^ (2.20) 

as follows from solving the {jip) Einstein equation (2J), using ( 2.19| ) to eliminate A4. 

More general solutions of type 3 cannot be found exactly, but an approximation is 
possible for values of E which are sufficiently close to the bottom of the potential in fig. 
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1(a) so that it can be treated like a harmonic oscillator. This gives 

M{p) ^ {zo + ecos ikp)f/^ (2.21) 

where zq is the position of the minimum of the potential, and k'^ = U"{zq). In general 
k and e are functions of Ag, A4 and (3. For example, in the case without magnetic flux, 
zq = (2A4/A6)^''^ and k"^ = ^Ag. If one wants to rescale the warp factor to unity, the 
physically meaningful place to do so is on the brane where the observer is supposed to be, 
for example zq + e ^ 1. 

Interestingly, in the more general solution ( ^.21 ), some dependence of the expansion 
rate on the tensions T3, Tg of the two 3-branes appears. An application of the jump 
conditions to be described in section 3.1 gives the following constraint: 

(2.22) 



rs - 27r [ zq- 



Tg — 27r \zq + 

This can be regarded as a fine tuning between the brane tensions required to obtain a given 
value of the expansion rate determined by A4. There is a one-parameter family of brane 
tensions which preserve the necessary relationship, showing that the expansion rate is not 
uniquely determined by the brane tensions. However the degeneracy of solutions along 
this curve in the Ts-Tg plane does not constitute a solution to the cosmological constant 
problem, and in this respect we disagree with the interpretation of In that reference, 
it was assumed that a deficit angle could appear spontaneously at the position of one of 
the branes, without having to specify any brane tension in the input to the construction 
of the model, namely the stress-energy tensor. Our point of view is that the 3-brane is 
not created by the geometry, but rather any singularity in the 2D curvature is due to the 
presence of a 3-brane with nonvanishing tension. (Another difference between our work 
and that of is that we do not assume the deficit angle, hence brane tension, to be zero 
at /9 = 0. This enlarges the class of solutions we are considering.) 

3. Brane tensions and jump conditions 

In the previous section we generalized the analysis of |^ to include extra sources of stress 
energy and to consider all possible signs of these various sources. A further new ingredi- 
ent which did not consider is the presence of branes which bound the bulk solutions 
given above. In each case we wish to allow for a 3-brane at the origin of the two extra 
dimensions, p = 0, which could provide support for the standard model and our observable 
4D spacetime. (A special case is a 3-brane with vanishing tension, which is the same as 
no brane as far as Einstein's equations are concerned.) The space will either terminate 
with another 3-brane or a 4-brane at some maximum value pm, or it may be noncompact, 
according to the possiblities mentioned in section 2, although the latter situation is not 
compatible with recovering 4D gravity at large distances. 

3.1 Jump condition at p = 



In all cases, the constant of integration R in (2^) is related to the tension T3 of the 3' 



brane at /? = 0. If is nonzero, there is a conical defect with deficit angle 5 given by 
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6 = KgTs, where Kg = SvrGg, in terms of the 6D Newton's constant. We can compute the 
ratio of the circumference and the radius of an infinitesimal circle around the origin as 
lime^o 27rL(e)/e = 27rL'(0). Therefore, continuing to work in units where kq = 1, we have 

T3 = 2^(1 - L'(0)) (3.1) 

If the space is terminated by another 3-brane at pm with tension fs, the relation becomes 

f3 = 27r(l + L'ipm)) (3.2) 

3.2 Metric jump conditions at p = pm 

On the other hand, if there is a 4-brane at pm, it may need to be accompanied by some 
additional form of stress energy in order to cut the bulk solution at this particular location, 
because there are separate jump conditions for the 69 and pp components of the Einstein 
equations. Therefore one generally needs two tunable parameters to satisfy both conditions. 
These conditions can be inferred from looking at the terms in the Einstein equations which 
have second derivatives or delta functions: 

M" IL" M'L' M'2 1 nc-/ ^ , . 

We will impose Z2 symmetry at the 4-brane so that the discontinuity in M' is —2M'. 
Integrating in the vicinity of the delta functions, we find 

ToO = 6^+2^; T/ = 8^ (3.5) 

where the functions are all evaluated at pm- If the stress energy tensor was given by a 
brane with pure tension, we would have Tq'^ = Tg^. This does not occur in the static AdS 
soliton solution ( |2.17| ), where M' /M — L' / L is always nonzero, except as p — > 00. One 
way to obtain a difference between Tq^ and Tg^ is to "smear" a 3-brane along the compact 
dimension of the 4-brane p3]. Since the 3-brane has Tg^ = 0, this procedure will give 



an extra contribution only to the T^^ components. The shift in Tq^ due to the smeared 

3- brane tension Tg is given by Tq" -T/ = T!^/[27rL{pm)]. This is not the only way to obtain 
a deviation Tq^ — Tq^\ for example the Casimir energy of a massless field confined to the 

4- brane gives — Tg^ ~ In general, the energy density of the source r' can scale 
like 1/L°^{pm)-, where a determines its equation of state. Conservation of stress-energy 
then dictates its contributions to the different components of the stress tensor [p[|,p^: 

T,^ = T, + V^{4>) + -^ (3.6) 

Tg' = T, + VM + (1 - o^)ji^^ (3.7) 

The scalar field potential Vm{(t>) is explained in the next subsection. 
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3.3 Jump conditions for (p 



In the absence of delta function sources for (/>, regularity of the solutions requires that 
4>'{0) = 4>'{pm) = 0. If potentials Vo{4')5^'^\p)^ Vm{4>)S{p — pm) are included in the La- 
grangian, (p' may be nonzero at these points. The effect of Vm is easily found by integrating 
the (p equation of motion in the vicinity of pm assuming Z2 symmetry: 

0'(.„.) = -i^| (3.8) 

The effect of Vq is more subtle. One can analyze it by regularizing the 2D delta function, 
S^^^p) — > l/(7rL'(0)e^) for p < e, and solving the field equation in the small p region 
using regular boundary conditions, (p'{0) = 0. The behavior of (p depends on the choice 
of potential. For Vq = \{4> — 4>o)'^, we find that (j) ~ ciJo{mp) + C2YQ{mp) in the interior 
region, with C2 ~ l/ln(e). Therefore in the limit e — > 0, we recover the boundary condition 
4>'{0) = despite the presence of the potential. On the other hand, the linear potential 
Vq = A(/> yields a singular solution, (j) ^ (j)Q ln(/3) near the origin. 



4. Stabilization in warped model 

Our ultimate goal is to understand the rate of expansion of the universe in models where 
we live on a codimension-two brane, with the expectation that the Hubble rate is due to 
the energy density of this 3-brane. Since the present universe has an energy density which 
is much less than that corresponding to the Planck scale, it is appropriate to treat A4 as a 
perturbation. We therefore want to first understand the model in the static limit, A4 = 0. 
Studies of 5D models have underscored the point that it is essential to stabilize massless 
volume moduli in order to recover 4D Einstein gravity at low energies, so we now consider 
a stabilization mechanism, in this section for the case where the 6D spacetime is warped 
by a negative bulk energy density. 

4.1 Stabilization by bulk scalar 

It has been shown [^2| that the radial size of the extra dimension is unstable in the AdS 
soliton model; one way to stabilize it is with a bulk scalar field which is prevented from 
vanishing by the presence of a potential 1^(0) on the 4-brane at p^. Following ||2^, we 
make the simple choice Vm{4>) = —^(p, which gives the boundary condition (p' = A/2, as 
noted in ( |3.8| ). In the large p region, where M ^ L g2fcp/5_ ^^^q large and small p 
regions, the solution for (p is approximately 

<^ = 0o(l + "iV/4 + ...), p-^l/k (4.1) 

with a± = —ki: and k = y^— SAg/S. The behavior at p = follows from 

assuming there is no potential for at p = 0. By matching the two solutions at p ~ 
and further assuming that m <C fc, one finds that B2 = 0{m^ /k^)(po, B\ = cpQ. Ignoring 
the B2 term, the jump condition at pm fixes a+(p{pm) = A/2. The assumption m <C /c is 
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useful for limiting the back reaction of the scalar on the geometry; in the limit m — > 0, 
(/)'—> 0, and the effect of the scalar is to shift the value of Ag, so that the solution for the 
metric would have the same form as that without the scalar. In this regime of parameters, 
the scalar solution can be approximated as (f) = (j)Qe^''^ ^/'^^ . 

In solving the jump conditions for the metric, it is useful to consider linear combinations 
which are independent of either the extra source r' or the bulk vacuum energy Ag. These 
linear combinations are 



X(t> + Ti = 






(- 


- 




a J 


ar' 


L' 


M' 


2L" 


T ~ 


1a 



= i-^i (4.3) 

where all quantities are to be evaluated at p = pm. The right hand side depends only on 
ratios where the overall normalization of the metric elements cancel out, whereas the left 
hand side depends on L itself. Using L = RM', V = RM", L'(0) = 1 - T3/2tt and ( ^ 
|2.17| ) to evaluate M"(0) = 2k'^/5, we can show that the constant R equals (5/2A;2)(l — ^). 
(Recall that we have rescaled the 4D coordinates so that M(0) = 1). Then L can be 
expressed as 

L = h-^ [l _ J) M' ^ k-^ (^1 _ ^) M (4.4) 



where the latter approximation holds for p ^ 1/k. Eqs. (^3, 4.3) determine pm, the 



position of the 4-brane, and they give a constraint on some function of T^, rs and r'. The 
latter is the fine-tuning which must be done in order to obtain a static solution. We will 
carry this out explicitly below. 

4.2 Back reaction of scalar (or magnetic flux) on metric 

To solve the jump conditions ( [4.2| , ^!^ ), we need to know how the scalar perturbs the metric. 
Although we will consider stabilization by the scalar and not the flux, for generality we 
indicate how the same analysis can be carried out for stabilization by flux. Similarly, we 
can also compute the perturbation away from the static solution by keeping the effect of 
the 4D cosmological constant, even though our immediate goal is to examine the static 
solution. 



Let 5p be the perturbation to M' /M in eq. (2.2), which when linearized takes the form 



V ^^,5, = \ (4^ - - m^^^ + ^ + ^) (4.5) 

The part of 6^ in the absence of flux, winding, or cosmological expansion, due only to the 
stabilization mechanism, is 



m 



2^2 fQ2{k+a+)p_i ^2a+p_i g2(~fc+a+)p 



+ 2—^ + 7—: (4.6) 



6Acosh^{kp) \ k + a+ a+ -k + a+ 
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It can be shown that in the m <^ k hmit in which we are interested, it is always true that 
0'^ <C m^cf^, so we dropped the derivative term to obtain ( 4.11| ). Similarly the decaying 



exponential part of (p gives only subdominant contributions in m^, so we approximated 
(j) = (poe'^+P. Now eq. ( [4 .61 ) can be integrated again to find the deviation in M itself, 

= f dpSfl^, ^ -^^E^f^Prri (at P = Prn) (4.7) 


where in the last approximation we have kept only the leading behavior in m?. The 
last result shows that \^kpm/m? must be kept small in order for the perturbation to be 
under control. Notice that we have used our freedom to impose the boundary condition 
(5M^(0) = to preserve the normalization M(0) = 1. 

To compute the analogous shift in L' /L, we perturb eq. (^^). In fact it will be 
convenient to compute Sl^f, — dfi^j,, 

6£^ -6p^ = -(a+ -) + -J-(</)'2 - mV') (4.8) 
\ pj tip 

Using £/p = (5 coth^(A;/3) — 3)/2 and p = (2/c/5) tanh(A;/j) for the unperturbed quantities, 
we can explicitly integrate ( [4.^ ) from to p. The result is 

This result will be essential for finding the size of the radial extra dimension p^ from the 
stabilization mechanism. 

4.3 Solution of jump conditions 

We are now ready to solve for pm and the fine-tuning between stress-energy components 
which leads to the static solution. Toward this end, it is useful to rewrite the difference 
between eqs. (|2.2D and (|2.l| ) in the form 

(M^L(;.-^))' = M^l(^ + !^ + ^) (4.10) 

If we are only interested in stabilization by the Goldberger-Wise mechanism, and the static 



solution, then the right hand side of (4.10) vanishes, leading to 



where we evaluated the equation at p = to evaluate the constant of integration c. Com- 
bining ( 4. lip with the jump condition (^^), we find the equation which determines the pm, 



namely 2L'(0)/ar' = M^L^~". The only hope for obtaining the desired large hierarchy 
Pm S> 1/k is in the case where a = 5. Otherwise we obtain g^^^""/^)'^''™- ~ L'(0)/r', which 
gives a large hierarchy only by tuning r' to be unnaturally small in units of the 6D gravity 
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scale. On the other hand, in the case a = 5, we must rely on the small correction (|4.9| ) to 
find the leading p dependence in the jump condition, 

The result for pm is 

k , / Am? , , 1 \ / . X 

In TTT(l-y) + n (4.13) 



m2 V5A2 ^ ' 2, 

This solution only exists if r' lies in the range 

y ~ ^\ <t'< ~ ^\ (4.14) 

This range could have had a large upper limit if SA^/Sm^ < 1. However we found in ( [4.7| ) 
that in fact we must have }? jm} < l/{kpm) in order for perturbation theory to not break 
down. So we see that there is still a fine tuning on r' for this model to give a large hierarchy. 
It is qualitatively different from the kind of tuning that would occur in models with a < 5; 
for these one would need to make the magnitude of r' exponentially smaller than its natural 
scale in terms of the 6D gravity scale. For a = 5, the magnitude of r' can be natural, but 
it must be very close to a specific value. One virtue of the scalar field is that it reduces 
the amount of tuning which is necessary. A similar tuning was present in the absence of 



the scalar field: (4.12) shows that y must be exponentially close to 1 if i;^ is absent, and if 
we insist on getting a large hierarchy. Once (j) is introduced, y must still be close to 1, but 
no longer with exponential precision. 

We could simply give up the tuning of r' and live with a model which does not try to 
solve the hierarchy problem. But since the main purpose of this paper is to explore the 
strange cosmological properties of this kind of model, questions of naturalness of the choice 
of parameters are secondary. It is convenient to adhere to the strongly warped case, since 
the hierarchy simplifies the algebraic analysis. We have seen that even in the completely 
unwarped case, similar behavior is revealed, concerning the unusual dependence of the 
Hubble rate on the brane tension is concerned. 

The remaining jump condition ( [4.2| ) is the fine tuning that must be imposed to make 
the effective 4D cosmological contant A4 vanish: 

k = br'e ^^P- 1 + --.kp^rn (4.15) 



2m? 5 \ IQm^ 

underscoring the fact that ultimately one needs an additional mechanism to fully solve the 
cosmological constant problem in these models. Our purpose here is not to address the 
problem of the 4-D cosmological constant, but rather the possible cosmological solutions 
of codimension-2 braneworlds. 

4.4 Stabilization by magnetic flux? 

We can attempt to stabilize the warped solution using magnetic flux; however since the 
effects of the flux are highly suppressed at large p because of the large exponential factor 
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it is not clear that it actually can stabilize the highly warped solution. Following 
the same steps as for the scalar field, we find 

^2 dpM^'^ = ^ ( r^^r. p-6fcp/5 _^ r^p-c^kp 

8 cosh {kp) Jq S/ccosh (kp) 



Scosh^C' ' 1-^^' V 

he approximate 
and in the large p limit 



where the approximate fit has cq = 1.38, ci = —1.91, C2 = — cq — ci, C3 = —(1 + |ci)/c2 



where do = —0.66, di = 5.5. The resulting equation for p^ has the same form as ( 4.12| ) 



except now the term 5Lp/L — 5Mp/M depends on pm only through g-^f^Pm^ which is of the 
same order as terms we ignored in ( [4.12| ). As explained there, this kind of dependence does 
not allow for a large hierarchy unless r' is tuned with exponential precision to a special 
value. This is exaclty the same situation as in the model with no stabilization, leading one 
to suspect that magnetic flux does not actually stabilize this model, unless the warping is 
weak (which would be the case if Ag > 0) . 

5. Priedmann equation in warped model 

We now come to the main point of the paper, in the context of the AdS-type models: the 
rate of expansion in the models we have described depends on the energy density (tension) 
of the standard model 3-brane in a strange way, relative to our expectations from 4D 
general relativity. Such an outcome would be less surprising if the low energy limit of the 
model was a scalar-tensor theory of gravity, in which the coupling of the extra massless 
scalar to the brane tension could explain the exotic effect. However we will note that the 
conclusion is valid even when moduli have been stabilized. We now demonstrate this in 
the warped case of the expanding AdS soliton model, (2.18), supplemented by a scalar field 



which stabilizes the radion by the Goldberger-Wise mechanism [36, 22 1. 

The main work has already been done in the previous section; there we found the back 
reaction to the metric by treating the scalar field as a perturbation. In the same way we 
can find the perturbation to the static solution due to a small rate of expansion by treating 
the A^/M'^ term to first order. We have a double expansion, in A4 and in (j)'^. We will work 
to first order in either quantity, but not including the mixed terms of order A^cp"^. 

The perturbation analogous to (|4.6D is 



5^,, = ^M-' l^p^M' (5.1) 

^ cosh-^(M {e'""^' - 1 - aie-^'' " 1)) (5.2) 

where a = (2^/^ — 1)^, b = (6/5) (2^/^ — 1)~^. The integral cannot be done analytically, but 
the appproximation (|5.2| ) is chosen to give the correct large p behavior, and it matches the 
exact value and its first derivative at /? = 0; it is accurate to 0.6% everywhere, as shown in 
fig. 3. 
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Figure 3: Ratio of tlie exact value of tlie integral in (4.11) to the approximation (5.2) as a function 
of kp. 



The shift in 6i is given by 



(5.3) 



By numerical evaluation and empirically fitting, we find that the following expression is a 
very good approximation to the integral of (|5.3D: 



(5.4) 



5L^ 6M^ 5A4 , , X 



This captures with great accuracy asymptotic approach of the solution to its large p limit, 
as well as having the correct small p behavior. The comparison between the exact result 
and the approximation is shown in fig. 4. 
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Figure 4: Numerical integration of (5.3) versus kp, and the approximation (5.4). 



The new contribution to the jump condition determining the combination i — fj, from 



( |4.10D gives 



^''°'^-5^^,M»-:,| ,5,5) 



where we used (4.4) to perform the integral exactly. This modifies ( [4.12| ) by giving a new 
term proportional to A4: 

where we have ignored all exponentially small corrections. 

Let us now recall what is the origin of the expansion of the universe according to our 
assumptions. We are interested in starting from a static configuration, and then perturbing 
it by adding some energy density on the 3-brane. Of course it is possible that the expansion 
of our own universe could be due to a small mismatch of bulk energy densities which are 
distributed in the extra dimensions. But our specific interest in this paper is to see how 
the energy density on the 3-brane affects the Hubble expansion. Therefore we will make 
the assumption that the quantity whose value is relaxed relative to the static fine-tuned 
situation is t^; hence in our perturbation series and A4 are of the same order. When 
comparing ( ^.6D to the corresponding relation ([4.12| ) for the static solution, there is one 
other quantity which must in general vary: the position of the 4-brane pm will deviate by 
some amount 6pm- However, all terms containing 6pm are or order A^cj)'^ or A| and so can 
be neglected. Therefore the difference between and ( [4.12| ) gives 



A4 = --^ —^^cosh "'{kpm) 

OTT 1 - 



2tt 



(5.7) 



where we used the constraint ( 4.14| ) to eliminate r'. This is the result we have been seeking, 
which shows an unexpected relation between A4 and 6T3: a decrease in the energy density 
on the 3-brane leads to cosmological expansion in this model. 

It is straightforward to give the more general version of eq. (^]^) , for arbitrary values 
of Q, and without assuming that we are perturbing around a static solution. We find that 

where M and L are to be evaluated at pm- This shows that the surprising dependence 
on T3 is not due to fine tuning of parameters. The Hubble rate is a decreasing function 
of r3 quite generally in this model. By choosing T3 so as to make A4 = 0, and perturbing 
T'i ^ T-},^ St3, one recovers eq. (5.7). 

When we perturb the other jump condition ( [4.2[ ), we obtain an equation for the shift 
in Pm, 

(^p' + m')) ^Pm + Sp^ + ^ (54 -6p^) = (5.9) 
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which gives 

_ -5A4 e'^P^ _ 3 5t^ 

32P coshl(A:p^)"20vrfer' ^'"'"^ 

6. Stability of unwarped model 



As aheady pointed out in eq. ( p. 20 ), it is possible to find unwarped solutions |13, 14 1 where 



the bulk is a two-sphere stabilized by magnetic flux. In these models the expansion rate is 
insensitive to the three brane tension. To be sure that this unusual behavior is not related 
to the presence of unstabilized moduli, in this section we perform a more complete analysis 



of the stability of the model than was done in ref. [30|. In the following we will set M = 1. 
We will consider in detail the static case and then give the main results of the analysis in 
the presence of the expansion. 

In the static case A4 = and, from eqs. ( |2.lD -( |2l^ ) we have 

^' + ^' = -/?', ^6 = Y. (6.1) 

The fluctuations of this model may be classified according to four-dimensional Lorentz 
transformations. On top of the usual transverse and traceless tensor (i.e. d^^hy = /i^ = 0) 
corresponding to five degrees of freedom there are three divergenceless vector modes (cor- 
responding to nine degrees of freedom) and seven scalar modes. Overall, the perturbed 
six-dimensional metric has 21 degrees of freedom. These 21 fluctuations of different spin 
transform under coordinate transformations parametrized by the infinitesimal shift = 
{e^,ep,eg). The infinitesimal shift along the four-dimensional spacetime can be decom- 
posed, in turn, as the derivative of a scalar and a divergenceless vector, i.e. = + Cfi- 
It is now easy to see that there are 3 scalar gauge functions (ep, eg and e) and one vector 
gauge-function. Of the seven scalar degrees of freedom, three of them can be gauged away 
by using the scalar gauge functions. Alternatively, the seven gauge-dependent scalars can 
be rearranged into four gauge-invariant scalar fluctuations. In the context of the study 



of six-dimensional Abelian vortices |40, 41] it was shown that a convenient gauge choice 
brings the line element in the form 

ds'^ = (1 - 2'4))r]^,ydx^'dx'' + (1 + 2i)dp'^ + L{pf{l - 12ip)d9'^ + 2TrL{p)dpde (6.2) 

where ijf^i, is the Minkowski metric. Since all the scalar gauge functions are completely fixed, 
no spurious gauge modes will appear. Furthermore, it can be shown that the scalar fluc- 
tuations appearing in (^^) obey the same evolution equations of the four gauge-invariant 



fluctuations which have been defined, in general terms, in |41|. 

The four scalar fluctuations of the geometry are coupled, through the perturbed Ein- 
stein equations 

6Rab = Stab, tab = Tab - ^tJ^Gab, (6.3) 

to the scalar fluctuations of the sources whose evolution can be obtained by perturbing to 
first order 

^A(^/\G\F^'')=0■ (6.4) 
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The two relevant fluctuations of the source are the fluctuations of the p and 9 components 
of the vector potential. These fluctuations will be donoted as Ap and Aq. The divergence- 
full part of the four-dimensional vector potential is also a scalar but it can be gauged away 
by using the ?7(1) gauge symmetry. 

Now that the gauge is fully fixed, we are ready to find the spectrum of fluctuations 
for the class of unwarped models defined by (6.1). In connection with the automatic 
adjustment of the solutions to cancel out any dependence on the tension of the 3-branes, 
we are specifically interested in the radial modes and, therefore, the possible dependence 
of the perturbed quantities upon the coordinate 6 will be ignored. Using ( |6^ ) the various 
components of eq. 



are 



fj. = u 
PP 

ee 

^ie 
ep 



3/? 



^T./) , 

2L ^ 



/5 



Aa 



vr' + 2£7r 
□vr = 0. 



(6.5) 
(6.6) 

(6.7) 

(6.8) 

(6.9) 

(6.10) 
(6.11) 



where □ denotes the four-dimensional D'Alembertian, —d1 + (9|. Eqs. ( |6.5[ )-(6.11) should 
be supplemented by the evolution equations describing the fluctuations of the sources which 
can be obtained by perturbing eq. (|6.4D to first order: 



(6.12) 
(6.13) 



J^'q - INq + /3L(4V^' + e' - + UAq = 0, 



DA, = 0, 



A' = 0. 



Notice that vr and Ap are decoupled from the other equations. The resulting system of 
coupled equations ( 6.11| ) and ( |6.13 ) implies that both vr and Ap are massless excitations. 
Furthermore, since Ap is constant, eq. ( |6.10 ) can be easily solved with the result that 



7i{p) = j^-l3Apj dpL{p). 



(6.14) 



However this mode is not normalizable, hence it is unphysical. In the action perturbed to 
second order the kinetic term for vr appears multiplied by only one factor of L. Thus the 
kinetic term for vr diverges as p ^ 0, in agreement with the findings of pl| . 

The remaining longitudinal degrees of freedom (i.e. •i/', ^ and Aq) satisfy a coupled 
set of linear differential equations whose eigenvalues determine the stability of the solution 
(|6.l| ). The boundary conditions to be imposed at the origin are [^ ]: 



(^(0) + e(o) = 0, v'(o) = o, ^;,(o) = o. 



(6.15) 
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If we now consider eq. (p. 51) we see that, for p ^ 0, eqs. (|6.15D imply 

</.(o) + v(o) = 



It is convenient to define two new variables: 

X -Y 



X + Y 



(6.16) 



(6.17) 



2 ' "2 

With these variables we will have that the linear combination 3 x (/i, v) + (0, 9) leads to 

X" + UX' + (□ - I3'^)X = (3^Y (6.18) 

while the difference of the {p, p) and {0, 6) components of the perturbed Einstein equations 
1 ^ ) and (^ gives 

X" - Y" - e{X' - Y') = DY (6.19) 



and eq. (6.9) determines the gauge field 



L 



Ae = -{X' + 2iY). 



(6.20) 



In terms of X and Y the boundary conditions ( 6.15| ) are now X'(0) = ^^'(0) and 1^(0) = 0. 

It is easy to verify that one solution to the coupled equations ( |6.18| , |6.19D matches the 
results of ref. pCl], with X = constant, Y = 0, and mass eigenvalue 



(6.21) 



which is the radion mass squared of this model. Although the analysis of [^] was in a 
4D effective theory, their result is exact since the solution has no dependence on the extra 
dimension. The constancy of the wave function shows that this must indeed be the ground 
state of the system ( |6.18 , 5.191 ), since nonconstant solutions should have higher energy and 
correspond to Kaluza-Klein excitations of the ground state. Nevertheless to be thorough 
we have conducted a numerical search for other eigenstates using the shooting method, 
and verified that ( |6.21 ) is indeed the ground state solution. 

The same analysis can be done in the case where the background solution has de Sitter 
expansion, as in eq. ( ^.20 ), rather than being static. The system ( 6.18 , 5.19 ) becomes 



X" + 3£X' + x(^D-/3^ + = + 12H^) Y 

X" - Y" - i{X' - Y') : 



4 



(6.22) 
(6.23) 



The form of the bulk solution is the same as when H = 0, but the mass eigenvalue 
generalizes to 



(6.24) 



showing that a sufficiently large rate of expansion destabilizes the compactification. 
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Let us now discuss the vector modes, which must be considered in a fully consistent 



decomposition of the metric and source fluctuations 6D. It was recently shown in [14| that a 
massless graviphoton field is present in the low-energy spectrum. In the following we restrict 
ourselves to the spectrum of the vector excitations with only radial dependence. The 
dependence of the perturbed quantities, considered in |4C], would presumably correspond 
to KK excitations of the lowest modes which we seek here. 

Denoting by and the divergenceless graviphoton fields appearing, respectively, 
in the p) and [n, 6) components of the perturbed metric, the equations for the coupled 
system of vector fluctuations reads: 



+ iV^ = 0, 
DV^ = 0, 

Z'' + £Z'+i'Z. + aZ. + 2pA' 



A';, + eA'^ + aA, 



P[Z'+£Z^]=0, 



(6.25) 
(6.26) 
(6.27) 
(6.28) 



where is the divergenceless fluctuation of the vector potential and eq. (6.28) follows 
from the perturbed component of the evolution equation of the gauge field. 

From eq. ( |6.26D it follows that is always massless in perturbation theory, as was 
argued in [^]. This conclusion, valid for radial excitations, follows from the cancellation 
of the contribution of Ag and the magnetic flux in the perturbed equations. However eq. 
( 6.25 ) implies that L^^. Like the tt mode in ( 6.14 ) this mode is not normalizable and 



therefore unphysical. 

To analyze the coupled system of excitations for and Afj_ notice that if the mode is 
massless, then eq. ( 6.2^ ) can be integrated once to flnd that A'^ = (5Z^. This allows us to 



rewrite the equation for Z^ in a decoupled form: 



.29) 



It is already known that in the case of the six-dimensional Abelian vortex there exists a 



normalizable vector zero mode [40, 41 1 (see also 
mode takes the form 



A^^I3 I Ldp', 



In our case, recalling (6T), the zero 

(6.30) 



Z^ ~ 



L 



which is the ground state solution of eqs. (6.27) and ( |6.28| ). This mode is normalizable 
since the canonical fields related to A^ and Z^ are simply the original ones multiplied by 
\/L. The localization of the gauge mode related to A^ was recently invoked [|^, 41, 42| 
as a mechanism in order to localize gauge fields in the presence of thick Abelian strings in 
six-dimensions. 

The masslessness of the graviphoton mode is not a concern, as would be a massless 
scalar mode, due to the fact that the graviphoton couples only to T^e, which vanishes on 
the 3-brane. Since the graviphoton does not couple to standard model matter, it is not 
phenomenologically constrained like a Brans-Dicke scalar. Moreover it does not seem a 
likely candidate for explaining the tuning of the bulk geometry in response to a 3-brane 
tension, since this is a deformation which looks like a purely scalar mode. 
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7. Discussion and conclusions 



We have made a comprehensive review of 6D cosmologies which contain at least one 3-brane 
where the standard model could presumably be localized. Emphasis was given to several 
solutions in which the rate of expansion of the universe is either insensitive to the tension of 
a 3-brane, eqs. ( 2.14 , 2.19 ), or else depends on ra in a way which is at odds with 4D general 
relativity eq. (5^). This is despite the fact that we are considering values of rs which are 
weh below the scale of compactification and the mass scale of moduli which could account 
for a departure from conventional gravity. If the theories under consideration contained a 
massless radion, for example, one would not be so surprised to see such behavior. However 
both of the models we have focused on do have stabilized moduli. The AdS soliton model 
has only a single radion among the fluctuations with azimuthal symmetry. It is massless 
for a = 5 and has an exponentially small negative (mass)^ for a < 5 in the absence of 
stabilization, but with the Goldberger-Wise mechanism it acquires a positive (mass)^ of 
order m'^k^'^ (j? {pm)e^'^^P"^ ~ ^2g-2fcpm/5 -pox realistic values of the hierarchy, this 

gives a radion mass of order MeV, which is safe because the couplings of the radion to 
the 3-brane are suppressed by the 4D Planck scale. And in the previous section we have 
confirmed the stability of the unwarped "football-shaped" model by doing an exhaustive 
study of its perturbation spectrum. 

The tensions of 3-branes in two 6D theories fail to contribute to the Hubble rate in 
the expected way at low energies, despite the fact that these theories have no additional 
massless scalars that might explain a deviation from 4D general relativity. Instead, the 
explanation of this apparent failure of decoupling lies in the special way in which a 2D 
manifold's internal geometry responds to the tension of a codimension two brane. This can 
be seen from a direct dimensional reduction of the 6D theory to four dimensions @ . 

Consider the metric (l^^) with a fluctuating 4D metric, gi^fiu{x^). When we integrate 
the 6D Einstein-Hilbert action over the coordinates yi of the internal space, we obtain 

A,eff = ^ vls^ j d^y Va^M^R^ - i^.es (7.i) 

V4,eS = y^i'yM^ ^V5i(-^^2+A6+ V + +^T3,.<5(2)(y-yi)^ (7.2) 

where i?„ is the n-dimensional curvature constructed from gn,ij,u', T3,i is the tension of the 
ith 3-brane, and . . . denotes contributions from other sources such as the bulk scalar field 
or a 4-brane. In the unwarped solution, where rs^i = T3^2 ^ ''"3) there are two cancellations 
which make R4 independent of T3. First, the curvature R2 has singular contributions that 
exactly cancel the delta function terms explictly involving the tensions. One can see the 
crucial difference between codimension one and two branes by examining their respective 
contributions to R2 p4|: 



R2 1 

— = 2t3 5^'^\y) H — T4 6^^\y — 2/4) -|- nonsingular terms (7.3) 

Kg 2 

The coefficient of the 2D delta function is such as to cancel the similar contribution to 
the action from the 3-brane's stress-energy, whereas this cancellation is not exact for the 
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4-brane. The second cancellation of the dependence on the 3-brane tension is this: rs 
enters implicitly through the volume of the extra dimension, ^/goJfi) = L'(0) = 1 — 73/2-11 
which multiplies the entire nonsingular part of the action, including the coefficient of i?4 
in (7.1), which gives the 4D Planck mass: 



Ml = \ I d\^2M\y) (7.4) 

These two factors of -L'(O) cancel out of the Friedmann equation, = V^^eff/SMp, so that 
the Hubble rate is determined completely by bulk quantities, independent of the brane 
tension. Notice that one reason we can make this argument so simply in the unwarped 
case is the fact that the bulk solution changes only through a rescaling of L as a result of 
changing the brane tension. 

When we try to make a similar argument to deduce the form of the Friedmann equation 
in the warped solution, it is not straightforward, because M and L do change from their 
static forms when the brane tension is changed, which induces a nonzero value of A4. It is 
still true that the delta function term r^S^'^^y) is cancelled by a singular term of exactly 
the same magnitude in ^/g2R2■ But now there are additional terms of order A4 from the 
perturbation to the bulk part of the solution. These did not appear in the unwarped case. 
The essential difference between the two models is that the jump conditions at the 4-brane 
bring in a new source of dependence on L, hence the 3-brane tension, which does not exist 
in the type 1 solution or in the unwarped type 3 solution. These observations were also 
made in p. 

The hope would be to solve the cosmological constant problem by using bulk super- 
symmetry to ensure A4 = 0, in a model where A4 is independent of the 3-brane tension. 
Then A4 will be insensitive to all quantum corrections to the vacuum energy which arise 
if the standard model is confined to the 3-brane. However to describe realistic cosmology, 
one needs to depart from pure de Sitter space and look for Friedmann-Robertson- Walker 
solutions. In fact, it is not obvious that solutions to the 6D Einstein equations exist when 
cosmological matter and radiation are added to the brane. Let us suppose the line element 
has the rather general form 

ds^ = -N'^df + M'^dx^ + B'^dp^ + L^dO'^ (7.5) 

where all the metric functions can depend on t as well as p, and consider the terms in the 
Einstein equations which can have delta function singularities: 

f M" M'L'\ L" ^,r,. ^ , ^ 

JM" M'L'\ IN" N'L'\ L" „,,,,,, , 

/Af" M'L'\ /N" N'L'\ , , 

In all the solutions we have considered, the energy density of the brane is given by e = 
~P = T3, and M' = N' = at p = 0, so that there is no singular part in M" or N". 
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To obtain e 7^ —p, it is necessary to have singular behavior not only in L" /L, but also in 
the other terms. However L"/L gives a 2D delta function 6^'^\p) ~ 6{p)/p by virtue of L 
vanishing like p. We can only get this kind of behavior from M" /M if M also vanishes, 
which is not physically sensible since the standard model requires a nonvanishing metric. 

The only way we see to admit a general equation of state is to give the 3-brane a finite 
radius po- A solution to the equation M" + M'/r = —eM5^'^\p) can be found if S^'^^p) is 
a regularized delta function: 

M( ) = i P< Po .^g^ 

I Joiapo) + apoJiiapo) InJ, p> po 



where a = y^T/JiTp^. For branes of codimension one, the interesting corrections to the 
Hubble rate come from terms like (M'/M)^, which are nonsingular and have an unam- 
biguous value proportional to e^. For the regularized codimension two brane, (M'/M)^ ~ 
(ep/pg)^ which goes like l/pg when averaged over the brane. Thus we expect to get results 
which are sensitive to the internal structure of the brane, unlike the clean predictions that 
came from the codimension one case. Of course, unambiguous predictions could be made 
in a specific theory, for example where the brane was a cosmic string defect coming from 
an Abelian Higgs model. Whether the interesting cosmological properties of codimension 
two branes will survive after regularizing them requires further investigation. 
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Appendix 

A. Perturbation to AdS soliton from small A4 

Let us write E = Eq + 6E, where Eq = az^ and 6E = —bzQ^^ is an energy shift relative to 
the static solution which allows us to keep the same starting value zq; then (E — U) can 
be written as 

E - U{z) = a{zl - z^) + h{z^'^ - zl'^) (A.l) 
and the h term can be consistently treated as a perturbation, even in the region p = 0, 



z = zq. The integral ( 2.11 ) gives 




5p{z)- (A.2) 



6/5 

° dz + 0{b^) 



1.368 + I {^f^ - )' + 0((f ^ » ^0 



(A.3) 



fc3 I 0.85 el/2 -0.18 e3/2 + 0.062 e5/2 + 0(e^/2)^ t = 1 + e 

where now we define k = \2a\^/^ . This can be inverted to find z(p) and M{p) to first order 
in b, 
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B. Cosmology of codimension one branes in six dimensions 



In section 2, we mentioned that the case where z'{0) ^ corresponds to having a 4-brane 
rather than a 3-brane at the origin, and that such a model is similar to RS in 5 dimensions. 



Here we will show this explicitly, following closely the formalism of |37]. 
We write the metric as 



-r?{r,t)dt^ + a^{r,t)dx'^ + h^{r,t)dr'^ + c^{r,t)de'^ 



(B.l) 



and include a scalar to stabilize the bulk. We will write the branes' stress energy tensors 
as: 



5(5r)diag(Vo + P*, 1^) -p*, Vb -p*, Vb -p*,0,Vb - pl) 
+5{h{r - R))dmg{Vi + p,Vi- p, Vi -p,Vi- p, 0, Vi - p% 



(B.2) 



where the perturbations p, p have arbitrary equation of state. We will expand the metric 
components around the static solution 



n(r,t) = 5(r)e-^i a{r,t) = a(r)ao(t)e-^i("'*) 

b{r,t) = 1 + Bi(r,t); c(r, t) = a(r)co(t)e-^i('^'*) 
(/)(r,i) = 00 (r) +(t)i{r,t) 



(B.3) 



where d and (f)Q correspond to the static solutions and the other terms are higher order in 
powers of p. 

B.l Boundary conditions 

The complications related to having a codimension-2 brane, and a 4-brane for which Tq ^ 
Tq are not present in this model, since we are simply dealing with a space cut off by two 
codimension-1 branes. We simply impose S1/Z2 symmetry, so that the boundary conditions 
will be 



a c 



a' d n' 
2- + - + - 

a c n 



3^ + ^ 
a n 



y=R 
y=R 

y=R 



2 
2 

b{R)K'^ 



\y=R 
y=R 
\y=R 



(B.4) 
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c' 












n' 
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a c 




n 




+ 


n' 


V a 




n 



y=0 
y=0 
y=0 



2 ^ ols/=o 
2 'ly=o 



6(0)^2 



\y=0 ■ 



(B.5) 
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We must also consider the boundary conditions for the scalar field equation, which are 

l(l)'\r=R = VI (B.6) 

B.2 Order 

The static solutions must solve the Einstein and scalar field equations expanded to zeroth 
order in p. These can be written as 

d 



-4? = ^HViM + U'i + A) + 6 f ? y (B.7) 
a 2 \ a J 



Using the superpotential method |3M|, one can find an exact solution 



a(r) 


= c(r) = e-^°(^) 


^o(r) 


= -kr H ^ ( 

5 16 V 











32 

The boundary conditions at lowest order in p impose: 

Vb = — ^ 1 - —K vqc 



^-2ekr _ ^ 



k \ 16 

V.= '-f(l- ^.^.gee--^ ) (B. 



where 



A = -|,. (B.10) 

It is clear that at the static level, this model is completely analogous to Randall-Sundrum. 
B.3 Order p 

We now look at the perturbed equations. We first define the following variables: 

^-1 = - (^3A[ + C[ + 4^1 ^) + «^0o</'i 

Ti = N[- A[ 

6 = - (sa; + N[+ AB,^^ + K^4>'ocf>i 

Ci = Oi - </'o0'i + '/'o ^1 (B.ll) 
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which appear naturally in the boundary conditions. It is simple to check that these are 
invariant under the following gauge transformations: 



y = y + f{y); /(o) = /(i) = o 

Ai^Ai- 'if- Ni^Ni- if 
a a 

Bi^Bi + f- 01 ^ + ct^'of 

Ci^Ci- if. (B.12) 



a 

These new variables obey the boundary conditions 

^2 ^2 



Tri|y=ij = y (P + P); ^i\y=Q = -^{p*+P*) 
ei|,=« = -y/; 6U = yf^ (B.13) 



Since we still haven't fixed a gauge, we will now choose to work in the stiff potential limit, 
i.e. (;/)! = [^]. This allows us to write the first order in p terms of the Einstein equations 
as : 



ao 



ao_\ ^ oo Co 



ao Co 



o2Ao 



oq \ ^ Oo Co _ Co 

ao j ao Co co 

2 



Oo _^ / Oo 
«o \ao 



3e 



2Ao 



Oo , / ao 



+ 

Oo V^o 



+ 



Oo Co _^ 1 Co 
ao Co 3 Co 



ao 



*i + 3(^ Ti + 

.«0/ 1/2 



(^1 - 6 



1/2 



The solutions are given by: 
^ri(r) = e^^"^'^) ^'|r=0 + 3 

Ti(r) = e5^«(^) T|,=o + 



T; - 5^[)Ti 
^1 - 5A[,6 



(B.14) 
(B.15) 
(B.16) 



-A'o iCi + 4^1 - 3Ti) + k'cP'o'Bi (B.17) 



(B.18) 



Oo \ _^ Oo Co 
Q-o J ao Co 



6(0 



„5Ao{r) 



C|r=0 + 3 



Oo _^ / Oo 

ao V«o 



Oo \ _^ Oo Co _ Co 

ao / ao Co co 



-3Ao(r)^^ 



-3Ao(r)^^ 



(B.19) 
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The Friedmann equations are then, using SttG 



2/(f e-3^oWdr' 



Oo 
Go 

Go 



+ 



■ \ 2 



SttG r 



3 

-4ttG 
SttG 



Oo Co 

Oo Co 



p* + p* + e 



ip + p) 



+ 



1 Oo Co _ 1 Co 

2 flQ Co 2 Co 



(p^ + e-5^o(iJy 



We can combine these to write: 

gOO Co _^ ^ _ ^'^^ 
ao Co 



— 3p^= + e' 



-5AoiR) 



{p-3p) + 2{pl + e 



Co 3 

We can then use the fourth Einstein equation to solve for Bi: 
1 



-5Ao{R)pe 



(B.20) 



(B.21) 



B 



p* - 3p* 



p. 



(F(r) - SA'^e^^^) + e 



-5Ao(R) 



p-3p-p^) F{r) 



where we have defined 
F = 



(B.22) 



(B.23) 



These results are very similar to those one gets in the Randall-Sundrum scenarion. 
There are however some differences. First, p and p denote (4 + l)-dimensional quantities, 
since the branes in this model correspond to 4-branes. If the compact dimension is small 
however, there will be no light Kaluza-Klein modes, and one expects to recover standard 
cosmology at low energies. 

Another difference is the appearance of terms involving the time variation of the com- 
pact dimension's scale factor. Such terms are observationally constrained to be quite small, 
so that they should not lead to appreciable effects on late-time cosmology once a stabiliza- 



tion mechanism for the size of the angular dimension is included [39|. 



References 

[1] P. Binetruy, C. DefTayet and D. Langlois, "Non-conventional cosmology from a 
brane-universe," Nucl. Phys. B 565, 269 (2000) [arXiv:hep-th/9905012]. 

[2] C. Csaki, M. Graesser, C. F. Kolda and J. Terning, "Cosmology of one extra dimension with 
localized gravity," Phys. Lett. B 462, 34 (1999) [arXiv:hep-ph/9906513]. 

[3] J. M. Cline, C. Grojean and G. Servant, "Cosmological expansion in the presence of extra 
dimensions," Phys. Rev. Lett. 83, 4245 (1999) [arXiv:hep-ph/9906523]. 

[4] V. A. Rubakov and M. E. Shaposhnikov, "Extra spacetime Dimensions: Towards A Solution 
To The Cosmological Constant Problem," Phys. Lett. B 125, 139 (1983). 

[5] C. Wetterich, "The Cosmological Constant And Noncompact Internal Spaces In 
Kaluza-Klein Theories," Nucl. Phys. B 255, 480 (1985). 



-28- 



[6] G.W. Gibbons and D.L. Wiltshire, "Spacetime as a membrane in higher dimensions", Nucl. 
Phys. B 287, 717 (1987). 

[7] S. Randjbar-Daemi and C. Wetterich, "Kaluza-Klein Solutions With Noncompact Internal 
Spaces," Phys. Lett. B 166, 65 (1986). 

[8] R. Sundrum, "Compactification for a three-brane universe," Phys. Rev. D 59, 085010 (1999) 
[arXiv:hep-ph/9807348] . 

[9] J. W. Chen, M. A. Luty and E. Ponton, "A critical cosmological constant from millimeter 
extra dimensions," JHEP 0009, 012 (2000) [arXiv:hep-th/0003067]. 

[10] T. Gherghetta and M. E. Shaposhnikov, "Localizing gravity on a string-like defect in six 
dimensions," Phys. Rev. Lett. 85, 240 (2000) [arXiv:hep-th/0004014]. 

[11] G. R. Dvali and G. Gabadadzc, "Gravity on a brane in infinite-volume extra space," Phys. 
Rev. D 63, 065007 (2001) [arXiv:hcp-th/0008054]. 

[12] O. Corradini, A. Iglesias and Z. Kakushadze, "Toward solving the cosmological constant 
problem?," arXiv:hep-th/0212101. 

[13] S. M. Carroll and M. M. Guica, "Sidestepping the cosmological constant with football-shaped 
extra dimensions," arXiv:hep-th/0302067. 

[14] I. Navarro, "Codimension two compactifications and the cosmological constant problem," 
arXiv:hep-th/0302129. 

[15] A. G. Cohen and D. B. Kaplan, "Solving the hierarchy problem with noncompact extra 
dimensions," Phys. Lett. B 470, 52 (1999) [arXiv:hep-th/9910132]. 

[16] R. Gregory, "Nonsingular global string compactifications," Phys. Rev. Lett. 84, 2564 (2000) 

[arXiv:hcp-th/9911015]. 

[17] Z. Chacko and A. E. Nelson, "A solution to the hierarchy problem with an infinitely large 
extra dimension and moduli stabilization," Phys. Rev. D 62, 085006 (2000) 
[arXiv:hep-th/9912186]. 

[18] H. Collins and B. Holdom, "The Randall-Sundrum scenario with an extra warped 
dimension," Phys. Rev. D 64, 064003 (2001) [arXiv:hep-ph/0103103]. 

[19] J. E. Kim, B. Kyae and H. M. Lee, "Localized gravity and mass hierarchy in D = 6 with 
Gauss-Bonnet term," Phys. Rev. D 64, 065011 (2001) [arXiv:hep-th/0104150]. 

[20] Z. Chacko, P. J. Fox, A. E. Nelson and N. Wciner, "Large extra dimensions from a small 
extra dimension," JHEP 0203, 001 (2002) [arXiv:hep-ph/0106343]. 

[21] P. Kanti, R. Madden and K. A. Olive, "A 6-D brane world model," Phys. Rev. D 64, 044021 
(2001) [arXiv:hep-th/0104177]. 

[22] C. P. Burgess, J. M. Cline, N. R. Constable and H. Firouzjahi, "Dynamical stability of 
six-dimensional warped brane-worlds," JHEP 0201, 014 (2002) [arXiv:hep-th/0112047]. 

[23] M. Giovannini, H. Meyer and M. E. Shaposhnikov, "Warped compactification on Abelian 
vortex in six dimensions," Nucl. Phys. B 619, 615 (2001) [arXiv:hep-th/0104118]. 

[24] I. I. Kogan, S. Mouslopoulos, A. Papazoglou and G. G. Ross, "Multigravity in six 

dimensions: Generating bounces with flat positive tension branes," Phys. Rev. D 64, 124014 
(2001) [arXiv:hep-th/0107086]. 



-29- 



[25] O. Corradini, A. Iglcsias. Z. Kakushadzc and P. Langfoldcr, "Gravity on a 3-brane in 6D 
bulk," Phys. Lett. B 521, 96 (2001) [arXiv:hcp-th/0108055]. 

[26] J. Louko and D. L. Wiltshire, "Brane worlds with bolts," JHEP 0202, 007 (2002) 
[arXiv:hep-th/0109099]. 

[27] D. K. Park and H. s. Kim, "Single 3-brane brane-world in six dimension," Nucl. Phys. B 650, 
114 (2003) [arXiv:hep-th/0206002]. 

[28] T. Multamaki and I. Vilja, "Warped and compact extra dimensions: 5D branes in 6D 

models," Phys. Lett. B 545, 389 (2002) [arXiv:hep-th/0207263]. "Cosmological solutions of 
braneworlds with warped and compact dimensions," Phys. Lett. B 559, 1 (2003) 
[arXiv:hep-th/0301168]. 

[29] E. Ponton and E. Poppitz, "Casimir energy and radius stabilization in five and six 
dimensional orbifolds," JHEP 0106, 019 (2001) [arXiv:hep-ph/0105021]. 

[30] S. M. Carroll, J. Geddes, M. B. Hoffman and R. M. Wald, "Classical stabilization of 

homogeneous extra dimensions," Phys. Rev. D 66, 024036 (2002) [arXiv:hep-th/01 10149]. 

[31] Y. Aghababaie, C. P. Burgess, S. L. Parameswaran and F. Quevedo, "SUSY breaking and 
moduli stabilization from fluxes in gauged 6D supergravity," JHEP 0303, 032 (2003) 

[arXiv:hcp-th/0212091]. 

[32] U. Gunther, P. Moniz and A. Zhuk, "Nonlinear multidimensional cosmological models with 
form fields: Stabilization of extra dimensions and the cosmological constant problem," 
arXiv:hep-th/0303023. 

[33] G. T. Horowitz and R. C. Myers, "The AdS/CFT correspondence and a new positive energy 
conjecture for general relativity," Phys. Rev. D 59, 026005 (1999) [arXiv:hep-th/9808079]. 

[34] F. Lcblond, R. C. Myers and D. J. Winters, "Consistency conditions for brane worlds in 
arbitrary dimensions," JHEP 0107, 031 (2001) [arXiv:hep-th/0106140]. 

[35] L. Randall and R. Sundrum, "A large mass hierarchy from a small extra dimension," Phys. 
Rev. Lett. 83, 3370 (1999) [hep-ph/9905221]; "An alternative to compactification," Phys. 
Rev. Lett. 83, 4690 (1999) [hep-th/9906064] 

[36] W. D. Goldberger and M. B. Wise, "Modulus stabihzation with bulk fields," Phys. Rev. Lett. 
83, 4922 (1999) [hep-ph/9907447]. 

[37] J. M. Cline and H. Firouzjahi, "5-dimensional warped cosmological solutions with radius 
stabihzation by a bulk scalar," Phys. Lett. B 495, 271 (2000) [hcp-th/0008185], J .M. Chne 
and J. Vinet, "Order Corrections to Randall- Sundrum I Cosmology," JHEP 02 (2002) 
042, [hep-th/0201041]. 

[38] O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, "Modehng the fifth dimension with 
scalars and gravity," Phys. Rev. D 62, 046008 (2000) [hep-th/9909134]. 

[39] J. M. Cline and J. Vinet, "Problems with time-varying extra dimensions or 'Cardassian 
expansion' as alternatives to dark energy," arXiv:hep-ph/0211284. 

[40] M. Giovannini, "Gauge field localization on Abelian vortices in six dimensions," Phys. Rev. 

D 66, 044016 (2002) [arXiv:hep-th/0205139]. 

[41] M. Giovannini, J. V. Le Be and S. Riedcrcr, "Zero modes of six-dimensional Abelian 
vortices," Class. Quant. Grav. 19, 3357 (2002) [arXiv:hep-th/0205222]. 



-30- 



[42] S. Randjbar-Daemi and M. Shaposhnikov, "A formalism to analyze the spoctriim of brane 
world scenarios," Nucl. Phys. B 645, 188 (2002) [arXiv:liep-th/0206016]; "QED from 
six-dimensional vortex and gauge anomalies," JHEP 0304, 016 (2003) 
[arXiv:hep-th/0303247]. 



-31 - 



